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ABSTRACT 



Context. Although the question of progenitor systems and detailed explosion mechanisms still remains a matter of discussion, it is 
commonly believed that Type la supernovae (SNe la) are production sites of large amounts of radioactive nuclei. Besides the fact that 
the gamma-ray emission due to radioactive decays is responsible for powering the light curves of SNe la, gamma rays themselves are 
fi of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the 

O I kinematics of these explosion events. 

I Aims. We study the evolution of gamma-ray line and continuum emission of SNe la with the objective to analyze the relevance of 

O observations in this energy range. We seek to investigate the chances for success of future MeV missions regarding their capabilities 

ji of constraining intrinsic properties and physical processes of SNe la. 

!/3 Methods. Focusing on two of the most broadly discussed SN la progenitor scenarios - a delayed detonation in a Chandrasekhar-mass 

, ^, white dwarf (WD) and a violent merger of two WDs - we use three-dimensional explosion models and perform radiative transfer 

simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of ^*Ni and have similar optical 

I properties that are in reasonable agreement with the recently observed supernova SN 201 Ife. We examine the gamma-ray spectra 

►^ with respect to their distinct features and draw connections to certain characteristics of the explosion models. Applying diagnostics. 



such as line and hardness ratios, the detection prospects for future gamma-ray missions with higher sensitivities in the MeV energy 
range are discussed. 

Results. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost 
direct connection of the emission of gamma rays to fundamental physical processes occuring in SNe la permits additional constraints 
concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV 
^T missions such as GRIPS will resolve all spectral details only for nearby SNe la, but hardness ratio and light curve measurements still 

^^ allow for a distinction of the two different models at 10 and 16Mpc for an exposure time of 10* s, respectively. The possibility to 

C^ detect the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN la detections in 

^~H the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays. 

^ Key words, supernovae: general - hydrodynamics - nuclear reactions, nucleosynthesis, abundances - radiative transfer - gamma 

rays: general - line: formation 



1 . Introduction wavelength approach extending to the gamma-ray regime (MeV 

energies) is desirable. 

While there is general agreement that SNe la are the result of The gamma-ray spectra of SNe la are dominated by the 
thermonuclear explosions of carbon-oxygen WDs, many ques- lines of the ^^Ni — > ^^Co — > ^^Fe decay chain which pow- 
tions concerning the progenitor and explosion scenarios still ers the observable display of SNe la. In contrast to studies in 
remain open (cf. Hillebrandt & Niemeyer 2000). This lack of the optical or infrared wavelength regime where the emissivity 
knowledge is contrasted by the relevance of SNe la for mea- is strongly dependent on the complex opacity structure which 
suring cosmological distances as well as their influence on the in turn depends on the atomic level populations and chemical 
chemical evolution of the universe and emphasizes the need of a composition, the emissivities in the gamma-ray regime are deter- 
more thorough understanding of the underlying physics. To con- mined from their branching ratios and radioactive half-lives and 
strain current explosion models as tightly as possible, a multi- a few rather simple interaction processes like pair-production. 
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Compton scattering and photo-electrical absorption. This makes 
the gamma-ray regime an ideal tool to obtain a more direct 
handle on the mass-velocity distribution of the explosion prod- 
uct (Milne et al. 2004). These very promising prospects have 
led to numerous theoretical efforts (mostly one-dimensional) 
in investigating the gamma-ray emission of SNe la in the past 
(see for instance Clayton et al. 1969; Clayton 1974; Ambwani 
& Sutherland 1988; Chan & Lingenfelter 1988, 1990, 1991; 
Burrows & The 1990; Muller et al. 1991; Hoflich et al. 1992; 
Kumagai & Nomoto 1997; Gomez-Gomar et al. 1998; Hoflich 
et al. 1998; Sim & Mazzali 2008; Isern et al. 2008; Kromer 
et al. 2010; Maeda et al. 2012). Although there have been first 
attempts of taking three-dimensional effects into account (e.g. 
Hoflich 2002), a fully three-dimensional treatment of the ex- 
plosion hydrodynamics as well as the radiation transfer calcu- 
lations has not been carried out in previous SN la gamma-ray 
emission studies. We determine, to our knowledge for the first 
time, the predicted gamma-ray emission for completely three- 
dimensional SN la explosion models. The three-dimensional ap- 
proach is important for a realistic description of the distribution 
of the radioactive isotopes and the surrounding ejecta material - 
both factors the gamma-ray emission of SNe la is sensitive for 
as shown in the following. Furthermore, a thorough investigation 
of Une-of-sight effects due to different viewing angles relies on 
a three-dimensional treatment of the explosion scenario. 

However, as of now the detection of SNe la at MeV ener- 
gies has not been practicable owing to the low sensitivities or 
non-existence of detection instruments in this energy range. In 
this paper, we explore the additional benefits from the analy- 
sis of gamma-ray spectra towards a more sound theoretical un- 
derstanding of SNe la and discuss the detection limits of pro- 
posed next generation gamma-ray observatories. For this aim to 
be achieved, we run full detector simulations of the proposed 
MeV satelHte GRIPS (Greiner et al. 2012). 

In this paper, we focus on two main branches of suggested 
SN la progenitor models: The explosion of a carbon-oxygen 
Chandrasekhar-mass WD and the super-Chandrasekhar-mass vi- 
olent merger of two carbon-oxygen WDs. While the first sce- 
nario marks the end of a WD close to the Chandrasekhar limit 
that accretes mass from a companion star through Roche-lobe 
overflow or strong stellar winds, the second is considered to be 
the result of two closely orbiting WDs losing energy due to the 
emission of gravitational waves and merging finally. The latter 
scenario has attracted a renewed interest in the last years since, in 
contrast to previous thoughts, the publications by Pakmor et al. 
(2010, 2011, 2012a) demonstrated that violent mergers of two 
carbon-oxygen WDs can directly lead to a thermonuclear explo- 
sion while the merger is still ongoing. As an extension to pre- 
ceding gamma-ray studies of SNe la, we discuss the gamma-ray 
emission of a violent merger for the first time. 

Both scenarios generate optical observables that are similar 
to those of normal SNe la (Mazzali et al. 2007; Kasen et al. 
2009; Blondin et al. 2011; Pakmor et al. 2012b; Ropke et al. 
2012), in spite of significant differences in the total mass and the 
ejecta structure. Further observable distinctions through SNe la 
gamma-ray spectra would therefore provide an additional handle 
on the progenitor channel question. 

The outline of the paper is as follows: In Sect. 2 we de- 
scribe the applied explosion models and the radiative transfer 
scheme. After a discussion of the resulting gamma-ray spectra 
and a study of the applicability of different diagnostic tools in 
Sect. 3, we investigate the chances of detecting gamma-ray line 
emission from SNe la in the near future in Sect. 4 and conclude 
with a short summary of the main findings. 



2. Simulation frameworl< 

2.1. Explosion models 

In the first simulation, we calculate the gamma-ray emission 
for the explosion of a near Chandrasekhar-mass (Mch) WD 
as a delayed detonation (e.g. Khokhlov 1991). After an initial 
subsonic deflagration phase of nuclear burning that produces 
mainly iron group elements, the explosion turns into a super- 
sonic detonation. The remaining fuel is mostly burned to in- 
termediate mass elements due to the prevailing lower densities 
caused by the energy release of the deflagration mode and the 
subsequent expansion of the star We use the NlOO model from 
the set of three-dimensional delayed detonation simulations car- 
ried out by Seitenzahl et al. (2013) with the thermonuclear su- 
pernova code Leafs. For a detailed description of the applied 
techniques we refer the reader to Reinecke et al. (1999), Ropke 
& Hillebrandt (2005), Schmidt et al. (2006), Ropke & Niemeyer 
(2007) and references therein. We chose the NlOO model since 
it produces optical observables that closely resemble those of 
"normal" SNe la (Ropke et al. 2012). The model is based on an 
isothermal non-rotating WD in hydrostatic equilibrium with a 
central density ofp^ = 2.9 x 10''gcm"^. During the explosion, 
10^ tracer particles reproducing the underlying density profile 
of the WD record the thermodynamic conditions. We then use 
the information provided by the tracer particles to calculate the 
detailed isotopic composition in a post-processing step with a 
reaction network of 384 nuclides (Travaglio et al. 2004; Ropke 
et al. 2006; Seitenzahl et al. 2010). The initial chemical composi- 



tion is assumed to be 47.5 ' 



12, 



C, 50.0% '"O and 2.5' 



"Neby 



mass, resulting in an electron fraction of Ye - 0.498864, which 
corresponds to a zero-age main sequence metallicity comparable 
to that of the Sun. With a kinetic energy of 1.45 x 10^'erg and 
a total mass of 1.40 Mq of the ejecta, NlOO produced 0.604 Mq 
of ^''Ni. Roughly half of the ^''Ni is located in the inner 0.3 Mq 
at velocities below 4, 000 km s"' , while the other half is more or 
less isotropically but inhomogeneously distributed within the re- 
mainder of the inner ~ 1 .2 M©, out to velocities of 12, 000 km s"^ . 
The second simulation is the violent merger of a 1.1 Mg and 
a 0.9 Mq WD of Pakmor et al. (2012b). Again, this model repro- 
duces the features of "normal" SNe la at optical wavelengths 
(Ropke et al. 2012). To model the inspiral and the merger, 
Pakmor et al. (2012a) used a modified version of the SPH code 
Gadget (Springel 2005). Both WDs are constructed from a total 
of 1.8 X 10^ particles of equal mass. After a relaxation phase, the 
distance between the WDs is slowly decreased according to the 
method of Dan et al. (201 1) until the first particle of the less mas- 
sive (secondary) WD reaches the Lagrangian point between the 
two objects. This triggers the actual simulation to start. During 
its progress, more and more material from the secondary WD is 
accreted and heated up on the surface of the primary, giving rise 
to the formation of hot spots and the ignition of carbon burn- 
ing. Following the guidelines of microscopic detonation simula- 
tions (Seitenzahl et al. 2009), a detonation is initiated in a hot 
spot that reaches a temperature of more than 2.5 x 10^ K and 
a density of about 2 x 10^ gem"-'. The mapping of the actual 
simulation to a uniform Cartesian grid with 768 x 768 x 768 
cells and a total box size of 4 x 10''cm is then used as initial 
state for a simulation of the detonation flame with the Leafs 
code, where the detonation is ignited at the cell with the high- 
est temperature. For more information about the simulation see 
Pakmor et al. (2012a,b). The detailed composition of the ma- 
terial is again calculated with the tracer particle method and a 
post-processing step using a large nuclear reaction network (see 
above). The merger model yields are based on the same initial 



chemical composition as the delayed-detonation model (47.5 % 
'^C, 50.0 % ""O and 2.5 % ^^Ne by mass). With an asymptotic 
kinetic explosion energy of 1.7 x 10^ 'erg, the merger model pro- 
duced 0.616 Mq of ^^Ni out of the combined initial mass of 
2.0 Mq. In contrast to the delayed-detonation model NlOO, the 
^^Ni is mainly found at velocities below ~ 10, 000 km s"' and it 
is much more asymmetrically distributed in the ejecta. This is 
mainly due to the delayed explosion of the secondary into the 
already burned remains of the primary, which excavates a region 
virtually free of iron group elements at low velocity (see figure 2 
of Pakmor et al. 2012a). 



2.2. Radiative transfer 

Using detailed abundance distributions obtained from the tracer 
particle method we map the explosion ejecta to a 50 x 50 x 50 
Cartesian grid and follow the emission, propagation and inter- 
action of the gamma-ray photons with the Monte Carlo radia- 
tive transfer code Artis (Sim 2007; Kromer & Sim 2009). The 
main principles of the radiative transport code can be summa- 
rized as follows: In contrast to Nature's way of quantization, 
the radiation field is divided into Monte Carlo quanta represent- 
ing indivisible parcels of energy, providing several advantages 
concerning the simulation technique (Lucy 1999). Initially, the 
quanta begin as so-called pellets of radioactive material rep- 
resenting ""^Ni, ^*Co, ^^Fe and "^'^Cr (other radionuclei are ne- 
glected in Artis since they are typically not important at early 
times). Upon decay, the pellets are converted to mono-chromatic 
gamma-ray packets with frequencies sampled randomly accord- 
ing to the respective branching ratios. These gamma-ray packets 
are emitted into randomly chosen directions under the assump- 
tion of isotropic emission in the comoving frame. Then, their 
propagation is followed in frequency, three-dimensional space, 
and time, until they leave the ejecta or are removed from the 
gamma-ray regime due to interaction processes. The basic in- 
teraction processes of gamma-ray photons are pair production, 
photoelectrical absorption, and Compton scattering, with the lat- 
ter being the most dominant in the encountered energy ranges 
(Milne et al. 2004). In accordance with the ratios of the cross 
sections of individual interaction processes to the total cross 
section, the occurrence of a certain interaction type is sampled 
randomly. Throughout the simulation, a positronium fraction of 
zero is assumed, meaning that positrons, e.g. from pair produc- 
tion or nuclear decays, annihilate in situ and directly lead to the 
production of two gamma-ray photons at 0.5 1 1 Me V (cf Milne 
et al. 2004). The escaping gamma-ray packets are binned in fre- 
quency, time and direction, contributing to the spectral evolu- 
tion of the gamma-ray emission from the explosion event. Light- 
travel time eff'ects are taken into consideration. A thorough de- 
scription of the employed Monte Carlo radiative transfer scheme 
and additional references can be found in Lucy (2005), Sim 
(2007), Sim & Mazzali (2008) and Kromer & Sim (2009). 



3. Results and discussion 

Even though the Chandrasekhar-mass delayed detonation and 
the violent merger model show distinct differences in total 
mass and ejecta structure, both models produce nearly the same 
amount of ^^Ni (roughly 0.6 Mq). In spite of their differences, a 
straightforward preference to one of the models cannot be given 
by comparing the simulated optical spectra to the measured for 
the test case of SN 201 Ife (Ropke et al. 2012). In this section, 
we investigate whether additional discriminating features can be 
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1.0 
Photon energy [MeV] 

Fig. 1. Spectral evolution of the gamma-ray emission from the delayed 
detonation (red) and the violent merger model (dashed black) for dif- 
ferent epochs after the explosion (indicated in the upper left of each 
panel). The spectral spread due to different viewing angles is shown for 
the maximum-light epochs of the two models in gamma rays (indicated 
in light red for the delayed detonation model in the third panel and in 
gray for the merger model in the fourth panel). Part of the effect is ob- 
scured by Monte Carlo noise in particular in the continuum. This prob- 
lem is largely removed by using hardness ratios and broadband light 
curves. The ^*Ni and ^^Co emission lines are indicated in the first panel, 
the flux bands Ci, Li and L2 discussed in Sect. 3.2 are depicted in the 
second panel. 



identified from a study of the gamma-ray emission of the two 
models. 



3.1. Gamma-ray spectra 

Fig. 1 shows the spectral evolution of the angle-averaged gam- 
ma-ray emission arising from the delayed detonation and the vi- 
olent merger model (photon fluxes are always normalized to a 
distance of 1 Mpc). For a discussion of the importance of view- 
ing angle effects, the reader is referred to the end of this section. 

The spectra are dominated by lines of the decay chain 
56^1 -^ 56(-Q ^ 56pg £,^g jQ jjjg different half-lives of ^^Ni 
(6.1 d) and ^^Co (11. 2d), ^''Ni lines which dominate at early 
epochs (e.g. at 0.812MeV or 1.562MeV) vanish in the spec- 
tra at later times, according to the decreasing amounts of ^^Ni. 
Then, the spectra are mainly formed by strong emission lines of 
^^Co and a continuum contribution due to Compton scattering 
of line photons. The optical depth to Compton scattering mainly 
depends on the column density of target electrons. Since the en- 
ergies of gamma-ray photons of ~ 1 Me V are much higher than 
the corresponding binding energies of electrons in atoms, nearly 
all electrons, regardless if bound or unbound, are accessible as 
targets of Compton scattering processes. As the ejecta expand 
with time, optical depths are reduced, leading to an enhancement 
of the lines with respect to the continuum. 

Besides these common properties of both explosion mod- 
els, there are also some pronounced differences. The gamma-ray 
emission in the delayed detonation scenario evolves faster than 
in the WD-WD merger, a fact that is also visible in a compar- 
ison of the bolometric gamma-ray light curves (energy range 
from 0.05 to 4.0 Me V) of the two models (see Fig. 2, upper 
panel). The delayed detonation produces a peak photon flux of 
1.82 X lO^^cm^s"' at 54.3 days, whereas the merger yields a 
maximum flux of 1 .43x 10"^ cm^ s"' at 75.7 days after the explo- 
sion. Since the masses of produced radioactive ^''Ni are nearly 
equal in both models, this is purely a consequence of the fact 
that the transport of gamma rays is, to first order, only sensi- 
tive to the column density of the material above the emission re- 
gion. I.e. the larger total mass of the WD-WD merger delays the 
gamma-ray emission and gives rise to a lower (angle-averaged) 
peak flux compared to the delayed detonation model. The con- 
vergence of the two light curves at late times in Fig. 2 reflects 
the similarity of the ''^Ni mass in both models. In the optically 
thin limit, the gamma-ray luminosity is given by 
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where fco is the Ufetime of ''''Co (Sim & MazzaU 2008). Thus, 
if the distance to the object is known, late-time measurements 
of gamma-ray luminosities can be used to unambiguously deter- 
mine the explosion yield of ^^Ni. 

Another distinguishing feature can be seen in the early-time 
spectra of the delayed detonation model at 20.1 and 34.9 days af- 
ter the explosion: While two prominent lines of ^^Ni can be iden- 
tified at 0.158 MeV and 0.270 MeV, these two lines are nearly 
totally degraded in the merger model and vanish in the back- 
ground of continuum emission. This effect can be explained 
by the energy sensitivity of the Compton cross section. Since 
the cross section decreases with increasing photon energy, espe- 
cially lines at low energies suffer from effective Compton down- 
scattering and the additional contamination of likewise down- 
scattered higher energy photons. Therefore, the two lines can 
only build up if a significant amount of ^''Ni is located at small 
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Fig. 2. Bolometric gamma-ray light curve (upper panel) and bolometric 
UVOIR light curve (lower panel) for the delayed detonation (red) and 
the violent merger model (dashed black). The light curve spread due to 
different viewing angles is indicated in light red and gray. The photon 
fluxes resp. the luminosities are normalized to a distance of 1 Mpc. 



optical depths (cf. Gomez -Gomar et al. 1998), a fact that directly 
connects the occurrence of low-energy ^^Ni lines to the distri- 
bution of the radioactive material. These different distributions 
of ^''Ni can be seen in figures 1 and 2 of Ropke et al. (2012), 
where it is shown that there is much more ^^Ni at higher veloc- 
ities in the delayed detonation model than in the merger model. 
Furthermore, there is less material surrounding the ^^Ni regions 
in the delayed detonation model and the corresponding column 
densities are therefore lower than in the merger scenario. This 
property is clearly mirrored in the evolution of the gamma-ray 
emission, but cannot be inferred easily from measurements in 
other wavelength ranges (see e.g. figure 3 in Ropke et al. 2012). 
The larger optical depths outside the ^''Ni region lead to more 
efficient Compton down-scattering and thus a softer spectrum of 
the merger model. This is further enhanced by the fact that iron- 
group elements are confined to lower velocities in the merger, 
leading to less photoelectric absorption than in the delayed det- 
onation. 

The influence of different viewing angles on the gamma-ray 
spectra is illustrated in Fig. 1 for the maximum-light epochs of 
the two models in gamma rays. While the strongest lines of 
the delayed detonation model do not show much variation, the 
asymmetric structure of the ejecta in the merger model leads to 
spectral features of varying magnitude, and for certain viewing 
angles a distinction from the spectra of the delayed detonation 
model can be very difficult. The low-energy range of the spectra 
is still the most suitable for distinguishing the two models: In 
the merger model, the higher column densities due to the larger 
total mass result in more effective Compton down-scattering of 
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Fig. 3. Peak line flux ratios of the gamma-ray emission from the delayed detonation (red) and the violent merger model (dashed black). The 
graph on the left depicts the line ratio R^ = F(0.847MeV)/F(2.598MeV) of two '*Co lines, the graph on the right illustrates the line ratio 
R2 = F(0.158 MeV)/F(0.847 MeV) of a ^'Ni and a '*'Co line. The flux ratios are normalized to the optically thin Umit. 
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Fig. 4. Hardness ratios of the gamma-ray emission from the delayed detonation (red) and the violent merger model (dashed black). The ratio 
Hi = Ci /Li is shown on the left, the ratio H2 = L[ jLi is displayed on the right. To illustrate the sensitivity to viewing-angle effects. Hi is shown 
for different lines of sight in the left panel. 



the gamma-ray photons in this energy range than in the delayed 
detonation model. This hinders the formation of prominent low- 
energy ^^Ni lines. The effects of observing the two explosion 
models at different lines of sight are also illustrated for the bolo- 
metric gamma-ray light curves in Fig. 2. As before, the asym- 
metric distribution of ejecta in the merger model results in a 
much larger spread of photon fluxes compared to the delayed 
detonation model. Nevertheless, the different times of maximum 
flux, as well as the differing early time evolution of the light 
curves, provide clear distinctive features that can be used to dis- 
criminate the underlying explosion models. These differences 
are again due to the different total masses involved in the two 
explosion scenarios. The larger mass of the merger model leads 
to more efficient photon trapping at earlier times, delays the rise 
of the photon flux, and results in a flux peak occurring at later 
times. 

As a comparison, we show the bolometric UVOIR light 
curves of the two models in the lower panel of Fig. 2. For many 
viewing angles, a distinction between the delayed detonation and 
the merger model is nearly impossible. This is contrasted by 
the behavior of the bolometric gamma-ray light curves: Here, 
in spite of viewing angle effects, especially early time measure- 
ments as well as the determination of the maximum fluxes are 
very promising. This again underlines the advantages of gamma- 
ray emission studies of SNe la. 



In the following, we discuss further possibilities to discrim- 
inate the two explosion models on the basis of angle-averaged 
spectra and present the fundamental distinguishing characteris- 
tics. As described in the previous paragraph, viewing-angle ef- 
fects can considerably complicate the process of drawing infer- 
ences about certain explosion scenarios, but distinctions are still 
possible. 



3.2. Line and hardness ratios 

In addition to the analysis of line fluxes and light curves, line 
and hardness ratios represent further diagnostic tools that can be 
utilized for the study of the gamma-ray emission of SNe la (cf. 
Hoflich et al. 1998; Gomez-Gomar et al. 1998; Sim & Mazzali 
2008). While the information inferred from absolute flux values 
is always restricted by how well the distance to the source is 
known, flux ratios are distance-independent and not subject to 
this limitation. In Fig. 3, two characteristic peak-intensity line 
ratios of the WD-WD merger and the delayed detonation model 
are illustrated. Following Sim & Mazzali (2008) we define Ri = 
F(0.847 MeV)/F(2.598 MeV) as the ratio of two ^^Co lines and 
R2 = F(0.I58MeV)/F(0.847MeV) as the ratio of a ^''Ni and a 
'^''Co line. 
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Flux ratios are affected by similar processes as discussed 
in the previous paragraph: Due to the energy sensitivity of 
Compton cross sections, line ratios are, at early times before the 
optically thin limit, dependent on the column density of electrons 
and therefore on the material above the region containing the ra- 
dioactive species. The ratio of two lines of the same radioactive 
isotope is therefore determined by the opacity ratios at the re- 
spective line energies. This is shown in the left plot of Fig. 3 
for the case of Ri. Since the radioactive material in the merger 
model is behind much more opacity, the line ratio shortly after 
the explosion is much lower. This is different in the delayed det- 
onation model: Here, more ^''Ni and hence more ^''Co at lower 
optical depths increases the line ratio significantly. It is thereby 
not crucial to choose two specific ^^Co lines for the examination 
of such an effect, but the extent of the diagnostic validity relies, 
due to the energy dependence of the Compton cross section, on 
a sufficient spread between the energy of the two selected lines 
(cf. Sim & Mazzali 2008). 

In contrast to Ri, R2 is a function that decreases with time. 
Since R2 is a ratio of two lines from different isotopes, this 
behavior reflects the difference between the half-lives of ^''Ni 
and ^''Co: While the ^^Co line at 0.847 MeV strengthens, the 
0.158 MeV line of ^^Ni fades away at later times, leading to 
smaller peak flux ratios. The information value of R2 is based 
on the following aspect: As discussed before, the ^''Ni present 
in the outer shells at small optical depths is the main source 
of the emerging 0.158 MeV line, whereas the emission of the 
0.847 MeV line originates from the total abundance of ^^Co in 
the ejected material, especially at later times. Therefore, a larger 
value of R2 indicates a larger deposit of ^^Ni in the outer lay- 
ers of the ejecta and the different distributions of this isotope in 
the two models are clearly reflected in the evolution of the line 
ratios shown in Fig. 3. As stated before (Gomez-Gomar et al. 
1998; Sim & Mazzali 2008), the relatively low energy of the 
0.158 MeV line makes it also very sensitive to photoelectric ab- 
sorption processes. Since photoelectric opacities are dependent 
on the compositional structure of the ejecta above the radioac- 
tive material, appropriate ratios of low-energy ^^Ni lines to ^*Co 
lines at higher energies are well suited to study the composition 
of SNe la. The success in finding distinct abundance features of 
course relies on the quality of the measured gamma-ray data and 
for this reason also on the distance to the explosion site (see also 
the remarks in the next section). 

A way to partially circumvent the need for highest quality 
gamma-ray data for detailed SNe la studies is provided by the 
application of hardness ratios. Instead of discrete line intensi- 
ties, the fluxes of broader energy bands are compared to each 
other. On the one hand, in this coarser method, some informa- 
tion is necessarily lost. But on the other hand, in contrast to other 
wavelength ranges, the relative simplicity of gamma-ray spectra 
and the rather small number of factors that influence gamma- 
ray emission processes make hardness ratios an important diag- 
nostic tool. Following Sim & Mazzali (2008), Fig. 4 shows two 
hardness ratios. Hi = Ci/Li and H2 - L1/L2, as examples for 
the delayed detonation and the merger model. C] denotes the 
energy band from 0. 1 to 0.3 MeV, where the main contribution 
comes from continuum emission by Compton down-scattering 
and photo-absorption processes. The energy bands Li (from 0.7 
to 2.0 MeV) and L2 (from 2.0 to 4.0 MeV) are characterized by 
pronounced lines of ^^Co, while the importance of continuum 
emission decreases with higher photon energies. 

Similar to line ratios, hardness ratios are mainly affected by 
the energy dependence of the Compton cross section. Since Hi 
directly mirrors the strength of the continuum emission relative 



to that of the lines, the violent merger produces higher Hi val- 
ues than the delayed detonation - a direct consequence of the 
radioactive material's location at higher optical depths. Similar 
to the line ratio R2, Hi is an indicator of different column den- 
sities of target electrons. Since photoabsorption processes can 
be significant in the energy band Ci, Hi also sei-ves to a certain 
degree as an indicator of different compositions. H2 represents 
the flux ratio of two line-dominated energy bands and is analo- 
gous to the line ratio Ri. It depends mainly on the ratio of the 
Compton cross sections in the two corresponding energy bands 
and is also suitable to discriminate our explosion models. Hence, 
hardness ratios offer an alternative opportunity to extract infor- 
mation from the gamma-ray emission of SNe la - particularly 
for more distant explosion events. 

The results in Fig. 3 and 4 show that line and hardness ra- 
tios can be used as additional diagnostic tools to distinguish the 
delayed detonation from the merger model on the basis of their 
gamma-ray emission. These ratios are directly linked to the dis- 
tribution of ^^Ni in the ejecta and, due to their dependence on 
the column density of target electrons above the radioactive ma- 
terial, they are sensitive to different compositions and masses of 
the outer ejecta layers. Even in light of possible viewing angle 
effects (see left panel of Fig. 4 as an example), hardness ratios 
have proven to be quite robust in distinguishing our two models. 
Together with broadband light curve measurements, they pro- 
vide the best chances for conclusions on certain model features 
especially for more distant explosion events. 

However, besides the changes in the gamma-ray emission 
due to different viewing angles, other variations within a certain 
explosion scenario (e.g. additional rotation and mixing effects, 
modifications of details in the explosion mechanism) can lead to 
a spread in gamma-ray fluxes that makes it certainly very diffi- 
cult to distinguish between different realizations of SN la explo- 
sions. But concerning our two very different models represent- 
ing two different progenitor classes, we are confident that the 
delayed detonation and the violent merger as described above 
will leave a unique imprint on the gamma-ray emission: In the 
case of the same ^^Ni mass in both models, a higher total ejecta 
mass surrounding the radioactive material in the applied violent 
merger is unavoidable. CoiTesponding to our results that we have 
discussed before, this has a systematic effect on the gamma-ray 
observables. As we will show in the next section, the prospects 
for discriminating our two models in case of a nearby SN la are 
very promising, but rely on the realization of new MeV satellite 
missions in the near future. 

4. Detection prospects 

Despite the importance of the MeV energy range for the study 
of different astrophysical objects and processes (cf. Greiner et al. 
2012; Summa et al. 201 1), the mission with the highest sensitiv- 
ities at these energies was the COMPTEL instrument aboard the 
CGRO satellite in the mid-nineties. One possible successor with 
an enhanced sensitivity of a factor 40 is the proposed GRIPS 
(Gamma-Ray Imaging, Polarimetry and Spectroscopy) mission 
(Greiner et al. 2009, 2012). The focus of this section is on the 
progress that could be made in the future concerning the de- 
tection of gamma-ray emission from SNe la, using the GRIPS 
instrument as an example. 

Being able to detect photons with energies from 0.2 MeV 
to 80 MeV, GRIPS is specifically designed for the search of 
gamma-ray bursts and blazars, the study of particle acceleration 
and radiation processes in a variety of cosmic sources, and the 
study of supernova explosion and nucleosynthesis mechanisms. 
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Fig. 5. Simulated detector response for the GRM instrument aboard 
GRIPS using an exposure time of 10^ s. Shown is a comparison of the 
simulated angle-averaged gamma-ray spectra of the delayed detonation 
model (red) and the merger model (dashed black) at 25 days after the 
explosion for different distances to the source (indicated in the upper 
right). 



The main instrument of GRIPS is the Gamma-Ray Monitor 
(GRM), a combined Compton scattering and pair creation tele- 
scope consisting of two separate detectors with an effective area 
of 195 cm^ and an energy resolution of 17 keV at 1 .8 MeV. The 
so-called tracker, made of silicon strips, is the first detector 
where the initial interaction of the incoming gamma-ray pho- 
tons takes place. Except for the entrance surface it is surrounded 
by the second detector, a calorimeter composed of LaBr3 scin- 
tillator material that allows for the energy determination of the 
secondary particles. In the case of a Compton scattering event 
in the tracker, the incident gamma-ray photon interacts with an 
electron, whose energy and position can be measured. The scat- 
tered photon is recorded in the calorimeter where its energy and 
interaction point can be reconstructed. With these data, it is pos- 
sible to calculate the direction and the energy of the incident 
photon. In the case of a pair creation event, the incident gamma- 
ray photon is converted into an electron-positron pair within the 
first detector. The original direction of the incident photon can 
be determined from the tracked directions of these two particles. 
The energy of the two secondaries, and therefore the energy of 
the incident photon, is measured with the help of both the tracker 
and the calorimeter 



GRIPS measurements of gamma-ray emission from SNe la 
were simulated using the MEGAlib software package (Zoglauer 
et al. 2006). This package uses the Geant4 toolkit (Agostinelli 
et al. 2003) to simulate the passage of gamma rays and other 
paiticles through the detector, taking into account possible par- 
ticle interactions and decays. The detector response to an inci- 
dent gamma photon is calculated from the energies deposited 
by the particles in sensitive elements (silicon strips and scin- 
tillator crystals), taking into consideration instrumental effects 
such as the energy resolution and noise suppression thresholds. 
MEGAlib includes a dedicated event reconstruction algorithm, 
which is applied to the events that pass the trigger criteria (co- 
incident hits in the silicon strip detector and the calorimeter). 
This algorithm either reconstructs the event as a valid Compton 
scattering or pair creation event, or rejects it. In the case of a 
Compton event, the incident photon direction is constrained to 
a circle around the direction of the scattered photon. To select 
the events compatible with the source position of the photons, 
we required the minimum distance between the circle and the 
source position to be less than 2°. In addition, if the recoil elec- 
tron left a track in the silicon detector, the Compton scattering 
plane could be reconstructed and its rotation angle around the 
direction of the scattered photon was required to be within 30° 
from the source position. The model of the GRM detector used 
in this study corresponds to the setup described in Greiner et al. 
(2012). 

Measurements of gamma-ray emission from most sources 
are strongly affected by the presence of radiation backgrounds. 
In our simulations we include the background from diffuse 
cosmic photons based on Gruber et al. (1999) and the back- 
ground from albedo photons produced in the Eaith's atmosphere 
(Mizuno et al. 2004; Ajello et al. 2008). In Boggs (2006) these 
two components were found to be the dominating background 
for a Compton telescope at a low-inclination low Earth orbit. The 
generated background photons were processed through the same 
simulation and selection procedure described above and then 
the reconstructed background events were added to the source 
events. 

In order to account for the specific sensitivity range of 
GRIPS, we define here the hardness ratio H* analogous to Hi 
as the quotient Cj/Li, where Cf denotes the energy band from 
0.2 to 0.4 MeV. In Fig. 5, the simulation results for the measured 
gamma-ray spectra are shown for the two models at 25 days after 
the explosion. The exposure time of 10^ s corresponds to roughly 
5 days in the all-sky scanning mode of GRM, a time interval that 
allows for reasonable studies of the time evolution. Accurate 
measurements of the line and hardness ratios become difficult 
for large source distances due to background fluctuations and 
limited event statistics in the high energy part of the spectrum. 
For the GRIPS mission, nearly the best discrimination between 
the two explosion models is provided by measurements of the 
H* hardness ratio. The simulation outcomes of the measurement 
of //j are depicted in Fig. 6. This ratio can be used to distinguish 
the models up to a source distance of 5 to 7 Mpc. Our simula- 
tion results for light curve measurements are shown in Fig. 7. 
For a supernova event at a distance of 1 Mpc, a very detailed and 
accurate light curve can be measured and constraints on certain 
explosion models are possible. According to our analysis of light 
curve measurements, GRIPS should be able to easily distinguish 
the two explosion scenarios for a source distance up to 10 Mpc. 

Using the gamma-ray spectrum of the delayed detonation 
model at 60 days after the explosion as input, we additionally 
carried out detector response simulations for different source 
distances and a longer exposure time of 10'' s (roughly 12 days 
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Fig. 6. Simulation results for hardness ratio measurements with the GRM instrument (exposure time 10' s) for our two explosion models. On the 
left, a source distance of 1 Mpc is assumed, the right figure shows the results for a source distance of 5 Mpc. Taking into account the sensitivity 
limits of GRIPS, C[ here denotes the energy band from 0.2 to 0.4 MeV. The solid lines show the results of ideal measurements without background 
fluctuations and statistical errors. 



800 



^ 600 

o 



400 



3 



200 




merger 

delayed detonation 



10 



o 
O 



+4-+ 



4+ 



+ 



HH-ii-, 



A 



,+' 



4. 



.+ 



'^"h 
+■ 



H 1 Hh" 



'^m ++ 



j^ 



merger 

delayed detonation 



20 40 60 80 100 

Time since explosion [days] 



20 40 60 80 

Time since explosion [days] 



100 



Fig. 7. Simulation results for light curve measurements (energy range from 0.4 to 4.0 MeV) with the GRM instrument (exposure time 10' s) for 
our two explosion models. On the left, a source distance of 1 Mpc is assumed, the right figure shows the results for a source distance of 10 Mpc. 
The solid lines show the results of ideal measurements without background fluctuations and statistical errors. 



in on-axis pointing mode). The detection capabilities of GRIPS 
concerning the gamma-ray emission of SNe la are illustrated in 
Fig. 8, where the simulated measurements after background sub- 
traction are depicted. While the strongest lines at medium and 
higher energies can still be resolved even for sources at 20 Mpc 
and more, the low energy range of the measured spectra from 
sources at larger distances is dominated by background fluctua- 
tions. Although GRIPS will provide a spectral resolution of all 
details in the energy range from 0.2 to 80 MeV only for SNe la 
at distances up to a few Mpc, the application of line and hard- 
ness ratios enables valuable spectral studies for even larger dis- 
tances. For an exposure time of 10^ s, our simulations show that 
the hardness ratio //j can be used to distinguish the two models 
up to 10 Mpc. With this longer exposure time, GRIPS light curve 
measurements can serve as distinctive marks of our models up 
to distances of 15 to 16 Mpc. Even explosion events in the Virgo 
cluster will be accessible for gamma-ray studies with such an 
instrument. 

The success of such investigations of course depends on 
the accuracy of the applied background models (see above). 
According to Ajello et al. (2008), an uncertainty of a factor of 
three seems to be quite reasonable. Although the chances for a 
clear discrimination between certain explosion models decrease 



considerably for larger source distances, the ability of GRIPS 
in proving the gamma-ray emission of SNe la at distances be- 
yond 20 Mpc will lead to a significant improvement of SNe la 
detection statistics in the gamma-ray energy range. Due to the 
low sensitivities of previous gamma-ray observatories, only one 
SN la has been detected up to now (cf. Milne et al. 2004). 

In cases of events comparable to SN 2011fe, the detection 
of gamma-ray emission should be easily possible with an in- 
strument like GRIPS. Our simulation results show that impor- 
tant constraints on three-dimensional explosion models can be 
drawn immediately. Since the limited sensitivity of GRIPS be- 
low 300 keV hinders a thorough examination of this energy 
range for sources at larger distances, the ^''Ni lines at 0.158 MeV 
and 0.270 MeV, which can serve as important diagnostic tools 
(see also section 3.1), are not so easily accessible. This prob- 
lem can be cured by combining the gamma-ray data of GRIPS 
with the data of the likewise proposed ASTRO-H mission (cf. 
Takahashi et al. 2010) that covers the respective energy range 
with sufficient sensitivity. A recent study concerning the detec- 
tion capabilities of ASTRO-H shows that distance limits com- 
parable to those of the GRIPS mission can be reached (Maeda 
et al. 2012). The realization of two missions with sensitivities 
and photon energy coverages like GRIPS and ASTRO-H in 



A. Summa et al.: Gamma-ray diagnostics of Type la supemovae 




0.5 1.0 2.0 3.0 

Photon energy [MeV] 



Fig. 8. Simulated detector response for the GRM instrument aboard 
GRIPS for an exposure time of 10'' s. Shown is the simulated angle- 
averaged gamma-ray spectrum of the delayed detonation model at 
60 days after the explosion for different source distances (indicated in 
the upper right). 



the near future would therefore offer unprecedented detection 
prospects for the gamma-ray emission of SNe la and equally 
support the efforts in reducing the parameter space of current 
three-dimensional model concepts. Especially the degeneracy 
concerning optical observables of different explosion models for 
the recent explosion event SN 201 Ife (cf. Ropke et al. 2012) 
could be certainly removed by a study of the corresponding 
gamma-ray emission as shown in the previous paragraphs. 

5. Conclusions 

Using three-dimensional simulations of a delayed detonation of 
a Chandrasekhar-mass WD and a violent merger of two WDs as 
test cases, we show that the calculated gamma-ray spectra are 
well-suited to draw conclusions on the different ejecta structure 
of the two models, but viewing angle effects always have to be 
taken into account. The direct link of gamma-ray spectra to the 
abundance of ""^Ni - the radioactive isotope that powers the radi- 
ation in all other wavelength bands - and their reduced complex- 
ity due to their straightforward connection to fundamental phys- 
ical processes make them a promising utility that can be used 
in a complementary way to other measurements. Being mainly 
sensitive to the column density of the material above the emis- 



sion region, the gamma-ray emission of SNe la is well suited to 
study the composition as well as the total mass of the ejecta. Our 
analysis of different diagnostic tools such as line and hardness 
ratios demonstrates that especially hardness ratios offer the best 
prospects to distinguish the two models at further distances. We 
also underline the value of following the evolution of gamma-ray 
emission over an extended period of time: Even low-resolution 
spectra of distant explosion events lead to characteristic light 
curve shapes that allow a discrimination of the two models. If 
the flux maximum can be obtained to a precision of about 5 days, 
bolometric measurements are sufficient to discriminate the two 
models. The prospects of success of course depend on the future 
development in the sector of detection instruments. Concerning 
the sensitivities in the hard X-ray as well as in the soft gamma- 
ray range, a major step forward will be done with the planned 
ASTRO-H satellite and the proposed GRIPS mission, allow- 
ing for observations of gamma-ray emission from SNe la up to 
20Mpc according to conservative estimates. Although a spectral 
resolution of all line features will only be provided for nearby 
explosion events, our simulated GRIPS observations show that 
hardness ratio and light curve measurements can discriminate 
our models up to source distances of 10 to 16Mpc, dependent 
on the applied exposure time. This enhances the number of tar- 
get candidates in the MeV energy range significantly and will 
allow for more meaningful gamma-ray emission studies of SNe 
la than up to now, including potentially unique opportunities for 
model discrimination. 
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